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ABSTRACT 


A Field Oriented Control (FOC) strategy based on a 24-Sector Vector Space 
Decomposition (24S VSD) technique used to control a Dual Three-Phase 
Permanent Magnet Synchronous Motor (DTP-PMSM) applied on electric 
ship propulsion prototype is presented in this paper. This machine is supplied 
by Dual Three-Phase Voltage Source Inverter (DTP-VSI). This study carried 
out on these multiphase machines has revealed that the large zero sequence 
harmonic current components on (z1, z2) subspace, constitutes major 
drawbacks; despite their advantages. The machine’s dynamic model is 
accomplished in three two-dimensional and orthogonal subspaces. In order to 
reduce more extra stator harmonic currents that produces losses, the current 
work is also concering of to consider the three different approaches for the 
sake of boosting the machine’s efficiency. Thus, the principle selection of the 
reference voltage vector and the time calculation method are presented in 
details. The proposed strategy effectiveness is validated by simulation 
results. 
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1. INTRODUCTION 

Nowadays, the multiphase machines provide special characteristics compared with the three-phase 
counterparts. These machines offer the power segmentation and ameliorates the reliability and availability of 
the drive system. In electrical machines, the winding factors can be augmented, and the torque ripple reduced 
by using multiphases. Therefore, one of the perspective structures of multiphase systems is a Dual Three 
Phase (DTP) motor fed by two three-phase voltage source inverters (DTP-VSI). 

The synchronous machines are widely used in industry, particularly the permanent magnet 
synchronous machines (PMSM). Thus, the DTP Machines notably PMSM drives has become very attractive 
in high performance control [1]. And DTP motor drives are widely used for electric propulsion ships, 
locomotive traction, electric/hybrid electric vehicles and aircrafts [2-6]. 
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Several Space Vector Pulse Width Modulation (SVPWM) techniques for six-phase synchronous and 
induction drive system have been studied and developed in [7-9], such as conventional SVPWM and vector 
space decomposition (VSD) theory. 

In order to control the voltage vectors components (z1, z2) subspace, a Field Oriented Control (FOC) 

strategy applying 24-Sector Vector Space Decomposition (24S VSD) is presented and established, which can 
be divided into 24-sector Modulation strategy. For six-phase induction machine, the conventional VSD 
technique was reported in [1, 10]. Then, the 12-sector Modulation is developed in [2]. 
According to the VSD control [1], the original machine’s characteristics model, is transformed into three 
two-dimensional orthogonal subspaces in the stationary reference plane, identified as the torque-component 
(a, B) subspace, the harmonic-component (Z1, Z2) subspace and the zero-sequence (01, 02) subspace. [11-15]. 
It should be noted, the fundamental components contribute to the conversion of the electromechanical 
energy. Whereas, the (z1, Z2) and (01, 02) current components do not contribute to the last conversion. 
Therefore, they are restricted by the stator resistance and leakage inductance [16-18]; also, the extra current 
harmonics only produce stator losses. To reduce these extra current harmonics, a Field Oriented Control 
strategy using a 24-Sector VSD (24SVSD) technique is proposed and simulated in this paper. Thus, three 
approaches have been also considered and developed. Finally, the simulation results show that the adopted 
approach reduces significally the harmonic components in (z1, Z2) subspace. 


2. MACHINE MODEL 

The structure of the DTP-PMSM windings fed by the VSI- DC sources, presenting spatially dual 
windings three-phase stator are shifted by 30 electrical degrees with two neutral independent points, as 
shown in three Figures 1, 2 and 3. 
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Figure 1. DTP-VSI feeding to DTP-PMSM Figure 2. Machine windings 





Figure 3. VSI-fed machine 
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According to the VSD theory [1], the transformation matrix can be expressed as follows: 


1 1 _1 V3 v3 0 
2 2 #2 2 
V3 V3 1 1 
CS =a o 
1 th V3 v3 
Poot = a (1) 
48 B 1 1 i 
2 2 Z 2 
1 1 1 0 0 0 
0 0 0 1 1 1 


Applying the transformation matrix T, the Dual Three-Phase motor can be represented in three two- 
dimensional decoupled mutually orthogonal planes (a-ß), (z1, z2) and (01, 02). 

According to the referenced assumptions in [19], the machine’s all components equations can be 
identified in stationary configuration by following relationship: 


[Vas] = [Rs]liae| + £ [Pag] (2) 

[Yee] = [Lag]lice] + Yom: [o> | 3) 

LV, Z1, = a [Rs] [iz] + [Lz] < [iz] (4) 

[z Z1, = [Lz] [izz] (5) 

[v 01, = [Rs] [ee + [Lo] £ [icis] (6) 
where: 

a atts Taia cos 20 oo maig sin 20 

m CENA sin 20 Tatta = tati cos 20 a 

where: 


- (La, Lq ) : Direct and indiret inductances; 

- (Lz , Lo ) : Self transformed conventional and mutual leakage inductances; 
- (Wp) : Flux produced by permanent magnet; 

- (0) : Rotor position angle. 


The VSD technique [7] allows expressed the machine components in three two submodels sach as: 
the order harmonics 12K*" + 1(K = 1,2,3,...) are located in the (a, B) subspace, the order harmonics 
6K + 1(K = 1,3,5,...) are transformed into the (zi, z2) subspace and the order harmonics 3K “(K = 
1,3,5,...) are transformed into the zero-sequence subspace (01, 02). 

The (a-B) current components contribute to the electromechanical conversion. However, the currents 
in (Z1, Z2) and (01, 02) do not contribute efficiently in the electromechanical conversion which produces 
stator losses. 

The transformation matrix must be used to represent the stationary (a, B) subspace in the rotating 
frame (d,q) which is shown in the following [20]: 


cos@ sin@ 0 
S “sin 0 cos | (l, : four-dimentional unit matrix) (8) 
0 0 I, 


The electric and mechanic equations for this model in (d-q) plane are expressed by: 


d- ieie A 
k|- [o mlel + v2 16" = 
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where: p is the number of pair poles. 


3. RESEARCH METHOD “FOC BASED ON VSD TECHNIQUE” 
3.1. DTP- SVPWM Strategy 

The Pulse Width Modulation (PWM) techniques have been the focus of many researches. The main 
goal of the PWM is to control the inverter output voltage and to reduce the harmonic content in the output 
voltage. To increase the output voltage, the space vector pulse width modulation (SVPWM) technique is 
applied. This technique is used for adjustable speed drives. By using sectors in SVPWM, the SVPWM can 
identify the position of reference vector and the switches can be operated by the mentioned sectors. Thus, the 
field orientation approach consists on need of the (a-B) subplane as being reference stator voltage vector 
produced by the control system. Then, the harmonics in (Zzı-z2) and (01, 02) subspaces produces only losses. 
So as to reduce these harmonics. The average voltage vectors produced in the two subspaces must be 
fixed at zero [21]. 

In the Figure 4, two Proportional Integral regulators are used in the current control structure which 
generates the direct and indirect voltage reference, especially Va and V,. Furthermore, the closed-loop 
speed regulator yields the indirect-axis current. The speed control can be achieved by controlling the q-axis 
current iq provided that the d-axis current component ig is kept at zero [22-24]. Therefore, the control 
strategy of id=0 is a relatively simple method in ordinary 3-phase PMSM control strategy. When the d-axis 
stator current winding id=0, the mechanic (11) becomes: 


Tom = Pig Ya = Piq Pom (12) 
Vie 
> 
I; Vi S(abc)1 ml 





N 
ml 
ja 24-sector Mr DTP- 
2 k = e SVPWM nee i rn PMSM | 
OL PC 


_ on eo be | BA 
I 
(khc)i 
la 
Coordinate (abc)2 
Transformation 
I q 


Figure 4. Machine’s field oriented control under VSD technique 


3.2. 24-Sector Vector Spae Decomposition (24SVSD) 
The FOC technique for the machine drive is shown in Figure 4. 
The relationship between phase voltages and device switch states is as follows: 


Va 2 oh 21 0 0 0 4f%a1 
Vbi =f 2 sq 0 0 0 {Sn 
Vo4 0 0 (13) 


Spe 


Ve2 Se2 


Where: - (S = Sa, Sbi, Sci, Sa2, So2, Sc2): Switching modes; - E: DC bus voltage. 
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By using the matrix relationship (13) and the transformation matrix T, voltages in each submodel 
are expressed as follows: 


Va Vat 
Vg Vbi 
Vz1 Vo 1 

=T 14 
V32 Vaz N ) 
Vor Vo2 
Vo2 Ve2 


Figure 3 shows the wingdings and DTP-VSI fed DTP-PMSM. A combinatorial analysis shows 64 
switching modes. Therefore, 64 different voltage vectors can be applied to the machine. Each vector is 
represented by a decimal number, corresponding to binary numbers, considered in the order [Sai Spi SciSa2 
Sb2 S.2]. Thus, the both figures 5 and 6 shows the schematic diagram of the space vectors based on VSD 
respectively in (a-ß) and (z\-z2) subspaces. 


B Null vectors at the origin: Z; 
07 56635 


AM 


Z1 





Figure 5. Diagram space vectors based on VSD in («a — p) Figure 6. Diagram space vectors based on VSD in (Z4 — Z2) 


Regarding the Figures 5 and 6, the (a-B) voltage vectors can be decomposed into four dodecagons 
with different magnitudes (from innermost to outermost: D1, D2, D3, D4) [25]. Their magnitudes are identified 
by the following relationship: 


128) 


Upi = J3 
1 
Un> = zE 
e (15) 
Up3 = BE 
[(2+V3) 
Ups =e 


The maximum amplitude (a-$) voltage vectors become the minimum amplitude (z)-z2) voltage 
vectors, whereas the others keep the same magnitude. 
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Figure 7. Diagram space vectors based on VSD having 12 maximum magnitude and 24 half-magnitude in (a-B) plane 


In Figure 7, the 12 voltage vectors and 24 ones having respectively maximum and half magnitude 
divide the (a-B) subspace into 24 sectors. The voltage vectors selected allows having the smallest amplitude 
vectors in the (z1-Z2) plane, which guarantees the minimum current harmonics content in the (z1-Z2) plane. 
Thus, leading to the decrease losses [3]. 

The combination of the maximum magnitude (a-B) vectors and the ones with half magnitude allows 
synthetizing the reference voltage vector Vyer (Figure 7). According to the location of the reference voltage 
vector Vef in the (a-B) subspace (Figure 5), the proposed strategy uses three vectors having the maximum 
magnitude and the fourth vector having half magnitude. In the sector Sı, the voltage vectors (37-36-52) 
combined with the fourth voltage vector (60) and zero voltage vectors (7-56) are employed. The time 
corresponding to each voltage vector is determinated by the following relationship: 


1 2 3 4 
Va Ve Va Velje Vefa 
Ve Ve Ve Velie V 
Va Va Va Va “1 0 
Ven Va Va Va 0 


where: 
e T;: 1s the sample time; 
e t,t ts and t; : are the times corresponding to the four voltage vectors; 
e VĚ: is the projection of k” voltage vector on the x-axis (x= a, B, Z1, Z2). 


The time interval attributed to zero voltage vectors is the remaining time: 
to = I; — (ti T t> + ta + t4) (17) 


In order to allow the occurrence of, at most, two transitions during each PWM period, three 
approaches, using the 24-sector strategy should be considered. In the first approach (approach (1) 
henceforth), the switching sequence concerns the position of the null vectors in the middle and at the ends of 
the switching sequence. In the second approach (approach (2) henceforth), the null vector is placed at the 
ends of the switching sequence. In the third approach (approach (3) henceforth), the null vectors is placed in 
the middle of the switching sequence (Figure 8). 
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Figure 8. Switching sequence in sector | according to the placement of the null voltage vectors 


(a) Approach (1). (b) Approach (2). (c) Approach (3) 


By applying the same reasoning to the other sectors, the times of four non-null vectors is given by 


the table below (Table 1). 


Table 1. 24-sector modulation vectors applying times (S1-S12) 
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Table 1. 24-sector modulation vectors applying times (S13-S24) 


SECTORS 
S13 Si4 S15 S16 S17 Sig Sio S20 S21 S22 S23 S24 
T ti -Tı T4 -T6 Ts -T10 Ts -Tj To -T12 T2 -T3 T7 
I t2 -T2 -Tı -T; -Ts T4 -Tio Ts -Thy Ts -Ti2 To -T3 
= ts -T3 -T5 Tı -To Te -T2 Tio = Ti T4 Tiz Ts 
S t4 -T4 -T3 -Tg Tı -T5 To -To Tio -T2 Th -T7 Ti2 


Projecting voltage vectors on the (a-f) and (z1-Z2)-axis in the relationship (16), the coefficients in 
Table 1 are given as follows: 


Ti 1 —V3 
T, v3 -1 V3 —1 
T3 0 2 
T, 2 — V3 =l 
5 1 —2+v3 
6 Tg —1 Vrefa 
T; B 2E a 1 PHA (18) 
Ig “pqs. 1-73 
Tg 1 2-V¥3 
Tio 2 0 
Tii V3 1 
Tiz 1 V3 


The relationship (17) makes it possible to determine the time allotted to zero voltage vectors. 
The Table 2 as follows gives the final switching sequences for all sectors. 


Table 2. 24-sector modulation switching sequences (S1-S12) 


SECTORS 
Si So S3 S4 S5 S6 S7 Ss So Sio Su Siz 
V Voi 7 7 0 0 56 56 63 63 7 7 0 0 
E Vi 37 39 36 4 52 48 54 62 22 23 18 2 
s V, 36 37 52 36 54 52 22 54 18 2 26 18 
O V3 52 36 54 52 22 54 18 22 26 18 27 26 
R V4 60 52 55 54 6 22 16 18 58 26 31 27 
S Voz 56 56 63 63 7 7 0 0 56 56 63 63 
Table 2. 24-sector modulation switching sequences (S13-S24) 
SECTORS 
Si3 Si4 Sis S16 Si7 Sis Sio S20 S21 S22 S23 S24 
V Voi 56 56 63 63 7 7 0 0 56 56 63 63 
E Vi 26 24 27 59 11 15 9 1 41 40 45 61 
r V2 27 26 11 27 9 11 41 9 45 41 37 45 
O V3 11 px 9 11 41 9 45 41 37 45 36 37 
R V4 3 11 8 9 57 41 47 45 5 37 32 36 
S Voz 7 7 0 0 56 56 63 63 7 7 0 0 


4. SIMULATION RESULTS AND DISCUSSION 

The simulation tests were accomplished by using a 3-KW machine applied on electric ship 
prototype, whose parameters are indicated in Table 3 [19]. 

In all the simulation tests, the load torque is kept at 15 N.m, the machine speed is fixed at 300 rpm 
and the inverter switching frequency is maintained at the same frequency 5 kHz. The figures 9 - 14 show the 
simulation results of the proposed strategy. 

As shown in Figure 9, in a period of 24 sectors, it is observed that 6 sectors remain with no 
switching for approaches (b) and (c). This allows a 25% reduction of the switching compared with approach 
(a) where there is a continuous switching mode. Thus, a significant reduction of the switching losses is 
accomplished. 
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Table 3. Machine’s main parameters 


Value/ unit 


Designation 
DC bus voltage: E 


400 V 


1Q 
8.5 mH 
8.5 mH 
0.175 Wb 
89.10° kgm? 
0.01 Nms/rad 


statoric resistance : R, 


d-axis inductance : Lg 


q-axis inductance : Lq 
permanent magnet flux: Ypm 


total inertia moment: J 
total viscous friction coefficient: f 


number of pair poles : p 
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Figure 9. Swiching command and corresponding sectors in terms of positioning zero voltage vectors 


(a) Approach (1). (b) Approach (2). (c) Approach (3) 


In Figure 10, the stator phase current resultant from the three presented approaches is extremely 
sinusoidal. The Figure 11 shows the harmonic current analysis significantly decreased to 1.05%. Thus, the 


harmonic analysis shows that the THD is slightly lower in approach (1) compared to the two other 
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approaches (2) and (3). It should be noted that these stator current components will only produce losses since 
they have no contribution in the elaboration of air gap flux. 
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Figure 11 Stator Current Harmonic Analysis according to 


positioning zero voltage vector 


Figure 10 Stator Phase Current according to positioning zero ce Dena cE OA DEO E A RTRaED 


voltage vector 
(a) Approach (1). (b) Approach (2). (c) Approach (3) 


To be noted that for the three approaches envisaged; THD current and switching number, the 
following table: 


Table 4. Swiching sequences evaluation 


Appropriate positioning of the zero 


Swiching Number THD current Fundamental (20 HZ) 
voltage vector 
Approach (1) 6 1,05 7,292 
Approach (2) 5 or 4 1,34 7,290 
Approach (3) 4 or 5 1,34 7,299 


Figure 12 shows that the stator currents in the subspace (z-z2) are too weak (amplitude 0.2 A) 
compared with the stator currents fundamental (amplitude 8 A). In Figure 13, the harmonic currents do not 
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affect the torque output performance. Figure 14 shows that the machine speed achieves its reference with 
satisfying static and dynamic performance. 
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(c) 
Figure 12. Stator Phase Current in subspace (z1, z2) according to positioning zero voltage vector 
(a) Approach (1). (b) Approach (2). (c) Approach (3) 
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Figure 13. Machine’s Torque Figure 14. Machine’s Speed 


5. CONCLUSION 

In this paper, a FOC strategy based on the 24-Sector Vector Space Decomposition (24S VSD) 
appropriate to the control machine drive system applied on electric ship prototype has been presented in order 
to reduce the current harmonics. The chosen four voltage vectors divide the (a-B) subplane into 24 sectors 
which is related to electromechanical energy conversion and each sector is m/12 rad. This modulation 
approach is required to efficiently eliminates the stator current harmonics in (z1-Z2) and (01, 02) subspaces. 
The simulation results show that the proposed technique allows the notable reduction current harmonics. 
Thus, it is confirmed that the suggested strategy offers lower harmonic current analysis. This explains by an 
effective reduction to 1.05%. Also, it is worth noting that a significant reduction in the stator current 
harmonics is associated with a 25% decrease in the switching losses corresponding to approaches (2) and (3), 
is observed at the level of the converter-machine combination. The results attest the effectiveness of the 
suitable modulation. Performing experimentation by an implementation on Digital Signal Processor (DSP) 
board of the proposed strategy will be the subject of the further works. 
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